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Abstract 

The confluent SUSY QM usually involves a second-order SUSY transformation where the two fac¬ 
torization energies converge to a single value. In order to achieve it, one generally needs to solve an 
indefinite integral, which limits the actual systems to which it can be applied. Nevertheless, not so long 
ago, an alternative method to achieve this transformation was developed through a Wronskian differ¬ 
ential formula [Phys Lett. A 3756 (2012) 692]. In the present work, we consider the fc-confluent SUSY 
transformation, where k factorization energies merge into a single value, and we develop a generalized 
Wronskian differential formula for this case. Furthermore, we explicitly work out general formulas for 
the third- and fourth-order cases and we present as examples the free particle and the single-gap Lame 
potentials. 
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1 Introduction 

1.1 Background 

The supersymmetric quantum mechanics (SUSY QM), factorization method or intertwining technique 
are equivalent algebraic methods used in quantum mechanics to generate new solvable potentials. This 
technique starts from known solutions (either particular or general) of an eigenvalue problem in quan¬ 
tum mechanics and develop a deformation of it such that the solution of a new problem is given in terms 
of the original solution plus some parameters from the transformation. 

In this way we are able to obtain parametric families of new potentials, from which usually one tries 
to find those without added singularities in order not to change the domain of definition and to ob¬ 
tain physically acceptable quantum systems, which otherwise would be just mathematically acceptable. 
After years of development, this method has become highly technical, with a lot of options such as: in¬ 
creasing the order (from first to kth), changing the type of transformation (real, confluent and complex), 
among others. For a review of this topic see Refs. (TJEI and references there in. 

In this work we deal with a specific case of higher-order SUSY transformations, the confluent case, in 
which several factorization energies converge to the same value ||3}0 . Taking such limit appropriately, 
this transformation leads to more flexibility for spectral design compared to the usual real case. Let us 
note that the standard method to perform this transformation requires to solve an indefinite integral, 
which is sometimes difficult to accomplish. Despite that, there has been a recent interest in this type of 
transformations, both in the study of its general properties |6}^ and in its applications ISl ITOtiTSl . 

In recent times an alternative differential method has been developed in terms of Wronskians, 
which allows us to perform a confluent transformation without the need of solving integrals and, there¬ 
fore, expanding the pool of potentials we can deal with. In the present work, we will see that this 
differential method is more general than it was previously realized. Here we will prove the general 
algorithm for the fcth-order confluent SUSY QM or fc-confluent SUSY QM [him (note that for fc = 3 it has 
been called hyperconfluent SUSY QM llU). Furthermore, we will explicitly develop two different ways 
to calculate both the 3- and the 4-confluent SUSY transformation. 
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Finally, we will apply these new formulas to two interesting systems. The first one is the free particle, 
whose 3-confluent SUSY transformation is already known and it can help us to check our results. The 
second is the single-gap Lame potential, whose higher-order confluent SUSY transformation has not 
been addressed before. 

1.2 fcth-order SUSY QM 

In the fcth-order SUSY QM El, one starts from the following intertwining relationship, 

H,B+ = B+Ho, (1) 

where Hq is the original Hamiltonian and H}- is the new one, still to be found. Usually Hq is exactly- 
solvable but in general whatever we know (or do not know) about this system will be inherited by the 
new one. Furthermore, the more we know about Hq, the more flexibility we have to perform a SUSY 
transformation. Both Hamiltonians define their respective potentials Vq, I 4 and Bj^ is a general fcth-order 
differential intertwining operator, all of them given by: 

Ho = -dl + Vo{x), (2a) 

H, = -dl + Vk{x), (2b) 

(2c) 

;=0 

where the bj{x) are functions to be found. Also and B^ = (B^)^ fulfill the following factorizations 


k 



(3a) 

k 

KBk=UiHk-ej), 

/=i 

(3b) 


with ej being the factorization energies for the transformation functions Uj. 

According to the usual SUSY QM, the new potential I 4 will be a deformation of the initial one Vq of 
the form ElEEZES: 

Vk{x) = Vq{x) - 2dl[]nW{ui,U2, ... (4) 

where the functions iij are k mathematical solutions of the original Schrodinger equation for Vq and W 
is the Wronskian. We use the adjective mathematical to refer to the solutions of the eigenvalue problem 
without the corresponding boundary conditions, which are the physical ones. In this work we will fix a 
specific domain only in the applications part, such that the theory will be as general as possible. 

From now on we will consider the special case of the SUSY transformation where all the factorization 
energies Cj in Eqs. (jS) converge to the same value, i.e., Cj = e for j = In this case, Eq. (|U 

still remains valid (see Theorem 1 in Ref. ||6|) but the Uj are now generalized eigenfunctions associated 
with the same eigenvalue e. In fact, usually obtaining a generalized eigenvector of higher rank is quite 
difficult but, as we will see, departing from a known eigenvector there is a simple way to solve all the 
generalized eigenvalue equations associated with the same eigenvalue ll4l lT9ll20l and to find the related 
Jordan block (also known as Jordan chain or Jordan cycle). 

Usually in a fc-SUSY transformation, up to k new levels could be added. These states could be math¬ 
ematical or physical according to their boundary conditions in the sense we have just described and only 
the ones that are physical are counted as added levels to the spectrum of the new system. Neverthe¬ 
less, due to the nature of the confluent case, where all the transformation energies converge to the same 
value, only one level could potentially be added, depending again if the state is physical or not. The 
new eigenstate is given by IIbIITbI: 
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W{UI,U2,. . 

W{ui,U2,...,Uk) 


(5) 


where fc > 2 in order to have at least two factorization energies converging and W{u) = u and again in 
the confluent case these functions Uj are generalized eigenfunctions. 

In the usual confluent SUSY transformations, this system is solved through indefinite integrals and 
from there a new family of potentials can be derived. In this work we will look for the general solution 
of Eq. (|D in a different way, through Wronskians and derivatives. 


1.3 Generalized eigenvectors and Jordan blocks 

The equation that defines a usual eigenvector ip can be written as 

{H-E)ip = 0, (6a) 

!/> 7^ 0, (6b) 

where H is any operator and E is its eigenvalue. 

There is also a generalized definition given by 

(H - E)’^ip = 0, (7a) 

{H-E)’^-^ip^0, (7b) 

where ip is known as a generalized eigenvector of rank k of the operator El and E is the generalized eigenvalue 
of H Hill. 

Starting from a standard eigenvector (of rank 1) we can form a chain of consecutive solutions for 

the generalized eigenvectors of rank k. In the simplest case, let !/^i be an eigenvector of rank 1 of H, 

i.e., {H — E)ipi = 0, then an eigenvector of second-rank must fulfill (Ef — E)^ip 2 = 0, which can be 
accomplished if we ask that {H — E)ip 2 = Therefore we get the following system of equations: 


(H-E)!/>i=0, (8 a) 

{H-E)^2 = ^i, (8b) 


which is known as a Jordan block, Jordan chain or Jordan cycle, since the matrix representation of H 
would be in this case 



i.e., a Jordan block. Generalized eigenvectors and their corresponding Jordan blocks of higher-rank can 
be obtained similarly and the matrix representation would he ak x k matrix with the eigenvalue E in 
the diagonal and 1 in the superdiagonal. Our final purpose here is to study the Jordan block of rank /c; 
in order to accomplish this, let us first study some low-order cycles and then address the fcth-order one. 

In the rest of the paper we will derive the differential Wronskian formulas for confluent SUSY trans¬ 
formations of orders 2, 3, 4 and k. Each of these cases will be developed in three steps: first, we will 
define the Jordan block of the corresponding length. Then, we will obtain the general solution of the 
generalized eigenvectors for the corresponding rank and from them we will finally obtain the differ¬ 
ential formula for the fc-confluent SUSY transformation. In Section we will apply these formulas to 
two examples: the free particle and the single-gap Lame potentials. We will finish this work with our 
Conclusions and an Appendix. 
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2 Differential formula for 2-confluent 

2.1 Jordan block of length 2 

Let us consider the following pair of generalized eigenfunctions of H, of first (mi) and second { 112 ) rank, 
associated with the eigenvalue e: 


(H-e)Mi=0 , (10a) 

{H — e)n2 = ui, (10b) 

defining a Jordan block of length 2. The functions Uj do not necessarily have a physical interpretation, 
i.e., they could be mathematical eigenfunctions of H. Eq. dlOab is a second-order homogeneous equation, 
therefore its basis has two linearly independent solutions. They can be taken as u and its orthogo¬ 
nal function v = u-^ defined by W(m, u) = 1 (this is for simplicity, indeed the only condition is that 
W{u,v) ^0). Without lost of generality, we can choose 

Ml = M, (11) 

and then 

{H-e)u = 0. (12) 

From W{u,v) = 1 we can obtain that 

We can easily prove that u(x) also solves the Schrodmger equation (fl^ by applying 3^ to Eq. ifTSll . i.e., 

iH-e)v = 0. (14) 


2.2 Eigenvectors of rank 2 

The second equation of the Jordan block dlObb is an inhomogeneous second-order differential equation 
and its general solution takes the following form 


“2 = U2 + 


,(P) 


(15) 


where is the general solution of the homogeneous equation and denotes a particular solution of 
the inhomogeneous one. We can use the same basis to describe the homogeneous solution of this second 
equation because we are dealing with generalized eigenvectors of the same eigenvalue IfT^ . Then, the 
homogeneous solution is given by a linear combination of our basis functions u, v, 


11 ^ 2 '^ = C\u + Div, 


(16) 


with Cl, Di G M. 

In order to find the particular solution let us suppose from now on that u and its parametric 
derivative with respect to e, d^u, are well defined continuous functions in a neighborhood of e. Hence, 
by deriving Eq. ifT^ with respect to e we obtain 


(H - e) deu = u, (17) 

where we assume that u = u{x, e), dedxU = dxdeU and H = H{x) ^ H{e). If we compare Eqs. dlObb and 
(II3, we can see that 

= deU, (18) 
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i.e., deU is a particular solution of the inhomogeneous equation we were looking for. Finally, the general 
solution of Eq. dlObb is given by 

U2 = Ciu + Div + d^u. (19) 

This is is the most general solution that closes the Jordan block of length 2. The apparent asymmetry 
between u and v due to the last term of this equation is produced by the initial choice of mi as m. Never¬ 
theless, there is no loss of generality since this is just a parameter choice that simplifies the equations by 
eliminating one term at the end (which for higher-orders will be quite useful). 


2.3 Formula for 2-confluent SUSY 

We can obtain a new Hamiltonian from Eq. (|4) by calculating the Wronskian of the two solutions iq, U 2 
of the Jordan block. This is straightforward, and the result is (it has 2! terms): 

W{uuU2) = Di + W{u,deU), ( 20 ) 


where we make use of W(m, v) = 1. Thus, the new potential V 2 {x) from the higher-order SUSY formula 
(|4} becomes 

14 (x) = V{x) - 2dl]n[Di + W {u,deu)], (21) 

which represents an alternative way to calculate V 2 {x) through the confluent second-order SUSY trans¬ 
formation. Furthermore, the new eigenstate is given by Eq. (|5]l as 


ip 2 ix,e) 


u 

Di -t W {u,deU) ' 


( 22 ) 


Note that these equations were first introduced by Bermudez et al. in Ref. ||5l, where the general 2- 
confluent case for the free particle and single-gap Lame potentials were also solved. Moreover, special 
cases for the free particle were previously addressed by Matveev ll2Tll and Stahlhofen 1221 without the 
constant Dj. We must remark that this constant is very important, as it characterizes a one-parametric 
family of SUSY partner potentials for each factorization energy e from which we can choose the ones 
without singularities, i.e., the meaningful quantum systems. Due to the Wronskian structure and our 
choice of mi = u, the constant Ci does not appear in the Wronskian formula in Eq. (l20l . therefore in the 
end we obtain a 2 -parametric family of potentials, with e, Di as parameters. 


3 Differential formula for 3-confluent 

3.1 Jordan block of length 3 

The Jordan block of length 3 is 


(H-e)Mi =0, 

(23a) 

{H-e)u2 = Ml, 

(23b) 

(H-e)M3 = M2. 

(23c) 


In order to solve it, we can follow an analogous procedure; thus, ui and U 2 are also those given by Eqs. 
(ITTl and (IT^ because they close the Jordan block of length 2, we only need to solve Eq. (I23cb . which we 
will do next. 


3.2 Eigenvectors of rank 3 


The eigenvector M 3 of rank 3 is also given as a linear combination of a homogeneous and an inhomoge¬ 
neous part 


M3 = M 


(ft) 

3 



(24) 
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The general solution of the homogeneous part can be expressed in the basis u, v because M 3 is also 
associated with the same eigenvalue e; therefore: 

= CiM + D 2 V, (25) 

(p) 

where C 2 , D 2 G IR are new constants. A particular solution can be found if we apply the operator 
to both sides of Eq. iflTll . leading to: 


(H-e)^ = ^eu. 

Then, replacing Eqs. dlTll and ll26t into Eq. (I23cll . we obtain: 

= CjdeM + DidgU + 

and therefore the general solution M 3 is given by 


U3 — C2U “h D2V “h ■ 


2 ' 


where Ci, Di are the same as in Eq. l[T9j. 


(26) 


(27) 


(28) 


3.3 Formula for 3-confluent SUSY 

The three functions mj, M 2 , M 3 that close the Jordan block of length 3 are given by Eqs. (ITTb . (fl^ and (|28] |. 
If we calculate their Wronskian we obtain (it should contain 3! terms): 


W(mi,M 2 ,M 3 ) =(CiDi — D 2 )'W{u,V,dell) + D\'\N {u, V ,d ^v) + D{\N{u,deU,deV) 

+ ^V^{u,v,dlu) + ^W{u,deU,dlu). (29) 

We can further simplify this equation by using the Wronskian identities that are found in the Appendix. 
In particular, we use Eq. ([76l to express the three-function Wronskians in terms of two-function ones, 
and we obtain 


W(Mi,Ii 2 , M3) 


D 2 - CiDi -h Di ]N{u,dev) -^N{u,dlu) 
— {Div + d£u)W{u,d£u) 


u — D^v — DidgM 


(30) 


Both of these equations can be useful in theoretical and numerical calculations. 
Einally we use Eq. (lUl with fc = 3 to obtain: 

V 3 (x) = Vo(x) -2a^ln[W(Mi,M2,M3)], 


(31) 


where either Eq. (l29l or Eq. ll30b should be employed to calculate a new potential from 3-confluent SUSY. 
Eqs. 1 I 29 I 1 and (l30b have 4 parameters, i.e., e, Ci, Di, D 2 . In this case C 2 is the constant that disappears due 
to the choice of mi . 

On the other hand, according to Eq. © the added eigenstate of the new potential V 3 (x) with eigen¬ 
value 6 is: 


ip 3 {x,e) 


Di -|- W {u,deU) 

W(mi,M2,M3) 


(32) 
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4 Differential formula for 4-confluent 


4.1 Jordan block of length 4 


Although the procedure is clear and we could already generalize it to fcth-order, we will also analyze 
explicitly the fourth-order method because we will use it in the applications in Section |6l The Jordan 
block of length 4 is 


(H-e)Mi = 0, 

(H — e)u2 = «!, 
{H-e)u3 = U2, 
(H-e)M4 = W3, 


(33a) 

(33b) 

(33c) 

(33d) 


where the solutions mi, W 2 , and U 3 are given by Eqs. (fTTIl , i fT^ and ll25t , while 14 will be calculated next. 

4.2 Eigenvectors of rank 4 

Once again, M 4 is composed by a homogeneous and an inhomogeneous part. The inhomogeneous part 
can be calculated from the parametric derivative of Eq. (I26l l, i.e.. 



(34) 


iv) 

and from here we obtain u^' and then the general solution M 4 as 



(35) 


where C 3 , D 3 G IR are new constants. 

4.3 Formula for 4-confluent SUSY 

Now we calculate the Wronskian of u-[, 112 , 1 / 3 , M 4 as given by Eqs. (fTTIl . IT^ . (l28ll and ll35ll , to obtain (it 
contains 4! terms): 

W(mi,M 2 /M 3 /M 4 ) =(CiDiD2 — D4C2 — D2 + D1D3) W(u, U, dgM, dgU) 




(36) 


Using again the Wronskian identities from the Appendix, in particular Eq. liTTb to reduce the order of 
the Wronskian, we get 


W(ui,U 2 /W 3 /M 4 ) = W(mi,M2)[W(mi,W4) +]N{u2,U2,)] - W^ (14,1(3), 


(37) 


7 






Wronskian differential formula for fc-confluent SUSY QM 


D. Bermudez 


where each of these terms is given by: 


W{ui,U 2 ) =Di + W{u,deu), (38a) 

W(mi, M4) + W{u2, M3) = C1D2 - C2D1 + D3 + C? W(m, d,u) + 3 CiDi W(m, a^u) + W(m, d^u) 

+ ^Di]N{u,d\v) + \]N{u,d\u) + D^W(v,d£v) + ^W(di;U,d^u), (38b) 

2 6 2 

W(mi,M 3 ) = D 2 + CiW(u,deu) + DiW(u,dev) + ^W(u,d^u). (38c) 

Moreover, from Eq. iSI with k = 4 we obtain 

V4{x) = Vo{x)-2dlln[W{ui,U2,U3,m)], (39) 

and using either Eq. ll36ll or Eqs. ll37t and (l38l we can calculate a new potential arising from 4-confluent 
SUSY. The constant C 3 does not appear in Eqs. (l36l and (l37l l, as we expected now considering previous 
cases. Therefore we have a 6 -parametric family of new potentials, with e, Q, C 2 , Di, D 2 , D 3 . 

The new potential V 4 (x) will have an extra eigenvector in e, which is given by: 


xpi{x,e) 


W(mi,M2,M3) 

W(mi,M2,M3,M4)' 


where W(mi, U 2 , M 3 ) is given either by Eq. ll29ll or Eq. ll30ll . 


(40) 


5 Differential formula for /c-confluent 


5.1 Jordan block of length k 

We can now foresee the structure of the solutions for a Jordan block of arbitrary length k. Next, we 
will prove our hypothesis using mathematical induction assuming our formulas are valid for k and then 
proving its validity for the k+4 case. 

Eirst, we need to obtain the k solutions Uj involved in Eq. (SJ. We start from a Jordan block of length 
k for the Hamiltonian H and the factorization energy e 

(H-e)Mi = 0, 

(H-e)u2 = Ml, 

{H-e)iik = u^- 

The Schrodinger equation for the seed solution m: 

(H-e)M = 0. 


(41a) 

(41b) 

(41c) 

(42) 


which is chosen from the general solution from Eq. (14lab as we saw previously. Then, let us obtain the 
parametric derivative of this equation with respect to e 


Deriving again 


(H — e)deU = u. 


/ u \ a 

(H-e)^ = deu, 


(43) 

(44) 


dlu 


and one more time 





Wronskian differential formula for fc-confluent SUSY QM 


D. Bermudez 


Now, we are ready to prove by induction the formula for the fcth derivative of u. Let us start from the 
induction hypothesis given by 

= (46) 

and apply on both sides to obtain the general formula for the index k + 1: 

(47) 

□ 


5.2 Eigenvectors of rank k 

Now we can solve the Jordan block of length k from Eq. (l4Tt . From the previous sections, we propose 
that the system is closed by = u and the following k — 1 functions Uj (/ = 2 ,..., k): 


Uj = ^j-i 
1=1 




/-I a/-l 

-ED,-- 

1=1 


+ 


d' 


i-l 


(/-!)! (;-!)!' 


(48) 


where, as usual, 3^/ = / and 0! = 1. Furthermore, this equation reduces to Eq. ifl^ for j = 2, the base 
case for our induction proof. 

Now, let us assume that this equation is valid for the fcth-order case and let us try to prove it for k + 1 
case. From Eq. l l48ll for k + 1, we have 


k 

“S:+l = ^k-l+1 

1=1 


(7^ 


k 


+ ^k-l+1 

1=1 


d[ 


M ■ 


(49) 


By applying the operator (H — e) to this frmction and using the Schrodinger equations for u and v given 
by Eqs. (fl^ and llT4t and the general formula for the parametric derivative in Eq. ll47t . we obtain 


k 


{H-e)uj,+^ = J^Ck-i+i 
1=2 


d[ ^u 

(7^ 


k 


+ ^k-l+1 

1=2 


(7^ 


d^-^u 


(50) 


where we make use of the fact that for I = 1 we obtain Eqs. llT2ll and (IT4t . Then changing the labels 
Z —)■ Z + 1 in Eq. l ISOll we have 


i.e.. 


{H-e)u^+i 


k-l 


JY ‘-k-l 
1=1 


d[ ^u 

(7^ 


k-l 


+ ^k-l 
1=1 


dg 

(7^ 


+ 


d^-iy 


{H-e)uk+i = Uk, 


which finishes our proof. 


(51) 

(52) 

□ 


5.3 Formula for fc-confluent SUSY 

Using the k solutions from Eq. ll48ll and inserting them into the Wronskian formula for SUSY QM of 
Eq. iSI we can obtain a Zc-confluent transformation to, in principle, any order k. These formulas become 
increasingly complicated; nevertheless, we must remember that we are performing a higher-order con¬ 
fluent transformation, in which k states converge to the same energy (Zc > 2), the corresponding integral 
equations for the usual method of confluent SUSY QM are more complicated to solve, and for several 
systems they cannot even be calculated. 
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Also, it is expected that the A:-function Wronskian that appears in this fcth-order transformation can 
be reduced as we did in the third- and fourth-order cases, because the basis fulfills the Schrodinger 
equations (fl^ and ifTit . which can be used to lower the order of the derivatives in two, i.e., 

d^^Uj dxUj, iij) (53) 

where f\ is a fimction to be found. Thus, the final Wronskian of k fimctions will only be dependent on 
the k functions themselves and their first derivatives with respect to x, as it was the case for the second-, 
third- and fourth-orders, shown in Eqs. ll20ll , ll30ll and ll37ll . respectively, i.e., 

W(Mi,M 2 ,...,Mfc) = f2{ui,U2,...,lik,dxUi,dxU2,...,dxUk), (54) 

where /2 is another fimction to be found. We must note that these Uj functions will ultimately depend 
on the parametric derivatives with respect to e of the seed solution u and its orthogonal function v up 
to {k — l)th-order, 

W(mi, U2, . . . ,Uk) = fziu,.. ., V,..., d^U, . . ., djV,dl~^dxv), (55) 

where is a different function. All the derivatives dx will be due two-function Wronskians in the 
variable x or 

W(mi,M 2 ,. . .,Mjt) = f 4 [W{d‘^{u,v},d^{u,v})], (56) 

where 0 < a, j6 < k — 1 are integer indexes, {u,v} means either one of the functions u, v and /4 is a 
function to be found. Furthermore, these Wronskians will define a {2k — 2)-parametric family of new 
potentials with parameters e, Ci,..., Q_ 2 / Di, ..., D^-^. These are the type of equations shown in Eqs. 
l|20ll, la and ||33. 

Furthermore, the new potential 14 (x) could have an extra eigenstate in e given by Eq. ||5]l. 


6 Applications 

In this section we will apply the third- and fourth-order formulas that we have obtained in this work 
to two interesting quantum systems. The first of them is the free particle, where both the differential 
and integral algorithms for the confluent SUSY transformation can be worked out. This will help us to 
test our ideas and compare our method with the previously known transformations. The second system 
is the single-gap Lame potential, for which the integral formulation is not easy to apply. As far as we 
know, this is the first application of 3- and 4-confluent SUSY to the Lame potential, in either integral 
or differential Wronskian formalism, because they involve elliptic functions, which are difficult to deal 
with. 


6.1 Free particle 


A free particle is not subject to any force so we can choose its potential as V = 0. In order to apply our 
algorithm we need two functions u, v such that W{u,v) ^ OVx G IR and for simplicity we will choose 
W{u,v) = 1. These functions represent a basis for the two-dimensional solution space associated with 
the equation {H — e)u = 0. We can choose 

(u,z;) = , (57) 


where e = —k^ and for calculating the derivatives we will use the chain rule as de = (d^x) d^- The 
domain of this problem is x G (—oo, oo), i.e., the full real-line. 

For the 3-confluent SUSY transformation, we take iq = u, we compute the other two generalized 
eigenvectors U 2 , M 3 and then calculate its Wronskian using either Eq. | |29] | or Eq.ll30t, leading to 


W(mi,M2,M3) 



8k^(CiDik 


D 2 K 


Dix)e“-e3“]. 


(58) 
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Figure 1: Three Poschl-Teller potentials and their corresponding added physical states at e obtained through 
3-confluent SUSY transformations. The plots correspond to e = —1,X2 = 4 (blue), e = —2,^2 = 0 (orange) 
and e = —3, X 2 = —4 (purple). (Color online). 


In addition, using Eq. lISTIl . the new potential is given by 


= - 


64fc3e2U(e2J:r _ Dik)[{l + Cifc 2 - fcx)Die 2 *^^ - 20 ^ 02 ^ - D 2 k?-e^^^ + 2 D 2 )c (-1 + Cifc^ - kx)] 


[e4“ - 40^x2 -h 8K2e2*^^(CiDiK - D 2 K - Dix)]2 


(59) 


From SUSY QMd, it is known that the 3-confluent transformation on the free particle potential where 
only one energy level is involved leads to the Poschl-Teller potential, but the expression in Eq. ll59ll 
does not look like the Poschl-Teller potential since it has three free parameters (Ci, Dj, D 2 ) instead of 
one. What happens is that we are performing a general 3-confluent transformation on the free particle 
and indeed we obtain more freedom to choose these parameters, however if we require that the new 
potentials should not have any singularities we have to restrict to Dj = 0 and consequently, this result 
does not depend on Ci. We have: 

,, 64D2K5e2“ 

^3(3^) f^ 2 Kx on „S^ 2 ' (^^) 


(e2K.r_8D2K3)2' 


and if we reparametrize D 2 in terms of X 2 as 


D 2 = — 


,-2kx2 


(61) 


then we obtain 

V 3 (x) = — 2 k 2 sech^[K:(x-h X 2 )], (62) 

which is indeed the Poschl-Teller potential with two free parameters: the factorization energy e = —k^ 
and the translation parameter X 2 . In Fig. [T]we can see three examples of potentials and their correspond¬ 
ing added physical states obtained from the 3-confluent SUSY. 

Therefore, even though we have obtained a more general family of potentials from the Uconfluent 
transformation, if we restrict the results to non-singular potentials we go back to simpler ones that are 
equivalent to those obtained by a first-order SUSY transformation. 

Now, in order to calculate the 4-confluent SUSY transformation we can use either Eq. ll36ll or Eq. (l37ll : 
as both calculations are easy in this case. After some algebra we obtain: 

W(wi,M2,M3, M 4 ) ^^ 4 ^ + ^ 2^1 “ CfDi - D 3 ) -h 8 k^x(CiDi — D 2 ) 

tP'x id / \ 

4k^(II) 2 — “h 2 .ID 1 KX — —h 2~2 — ^2 — ]D\X^j 

DiX ]j3q-2kx 

+ ^(CiDi - D 2 ) + C 1 D 1 D 2 + D 1 D 3 - C 2 D 2 - d2 + . (63) 
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From here, we can calculate the new potential V 4 (x); the result is prohibitively long to be written 
here, but it can be easily calculated using a symbolic computation software. Nevertheless, even though 
we obtained a 6 -parametric family of potentials only a subset of them is non-singular. Furthermore, all 
of the non-singular sub-families we were able to find can be rewritten as Poschl-Teller potentials. 

In that case, the Poschl-Teller potential should be included in this more general family. In order to 
prove this, let us take the Wronskian from Eq. (l63ll and simplify it by taking the sub-family of potentials 


with Cl = Di = D 2 = 0 , then 


- (e2>cx_32D3K)2' 


(64) 


and similarly to the 3-confluent case, we reparametrize D 3 in terms of X 3 as: 


q-2kx3 


(65) 


thus we can obtain 

V 4 (x) = —sech^[K:(x -|- X3)], ( 66 ) 

which is again the Poschl-Teller potential. 

It has been proven that if we apply the 2-confluent SUSY transformation to the free particle we obtain 
the Poschl-Teller potential if we restrict ourselves to the non-singular case (for the differential algorithm 
see Ref. 0 , for the integral one see Ref. 0). Here we used the differential algorithm we have just 
developed to calculate the 3- and 4-confluent SUSY transformations from the free particle and we also 
obtain Poschl-Teller potentials for the non-singular case with both methods (see Fig.[l]l. It is known that 
a non-singular SUSY transformation for the free particle that adds a single new level to this system will 
always lead to a Poschl-Teller potential, a one-soliton KdV potential Il4ll5l l23ll24ll , and this is the case for 
the k-confluent transformation, so our result agrees with the previous ones. Nevertheless, it does not 
have to be this way for all potentials, in fact, the outcome for the single-gap Lame potential in the next 
Section points out to a different direction. 


6.2 Single-gap Lame potential 

The Lame periodic potentials are given by l25lf27ll : 


V{x) = n(n + l)m sn^(x|m) = n{n 

where sn(x|m) is a Jacobi elliptic function and p(x) is 
single-valued doubly-periodic functions I281I29II with 


1 ) 


p{x + cv') -h -{m + 1) 


(67) 


the Weierstrass elliptic function, both of them are 



d0 

(1 — m sin^ 0 ) 1 / 2 ' 


( 68 ) 


and co' = iK{l — m) the real and imaginary quarter-periods of the former and half-periods of the latter. 
The potentials defined by Eq. (l67l have 2n -|-1 band edges which define n + 1 allowed and n + 1 forbid¬ 
den bands. The Lame potentials have been used to model a non-relativistic electron in periodic electric 
and magnetic field configurations which produce a ID crystal Il30ll . Furthermore, they are particular 
cases of the associated Lame potentials that have been studied in the context of higher-order SUSY QM 
by Fernandez et al. 1311 . In addition, these potentials are described by an exotic SUSY structure with four 
supercharges and two bosonic generators I231l32j[35l and they also display hidden symmetries related 
to the soli ton d 5 mamics II36II37I . 

In this work we shall deal with the simplest case for n = 1, known as the single-gap Lame potential. 
The spectrum of the associated Hamiltonian is given by: 


Sp(H) = [m, 1] U [1 -h m, 00 ), 


(69) 
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Figure 2: Spectrum of the Lame potential with n = 1. The white bands correspond to the allowed energy 
region, a semi-infinite [1 + m, oo) and a finite one [m,V\. The dark region corresponds to the energy gaps, a 
semi-infinite {—oo,m) and a finite one (1,1 -|- m). 


i.e., it is composed of one finite energy band [m, 1 ] and one semi-infinite [1 -|- m, oo) (see the white region 
in Fig. |2]l. The structure of the resolvent set of H is similar, namely, there is a semi-infinite forbidden 
energy band (—oo,m) plus a finite one (1,1 + m) (observe the dark zone in Fig. |2) . In solid-state physics, 
this system is very useful to describe the electron dynamics in semi-conductors. The lower allowed 
energy band is called valence band and the upper one is the conduction band-, also the finite forbidden 
energy band is called the band gap (see Ch. 5 of Ref. Il38l l. 

To implement the confluent second-order SUSY transformation, the appropriate seed solution ui 
has to be associated to a factorization energy e\ which is inside one of the forbidden energy gaps, i.e., 
the dark zone in Fig. |2l and such that W(mi,M 2 ) 7 ^ 0 V x G R all this following the confluent SUSY 
formalism. Note that in this case, the domain of the eigenvalue problem is also x G (— 00 , 00 ). For our 
example this can be achieved by choosing ui as one of the two Bloch functions associated to e I261I27I . 
i.e., 

^ a[x + co' + 5) 

a{x + w') 

1 £r(x + a;'-^) 

^ o-'^{5)dsp{S) a{x + Lo') ^ 

where a and ^ are the non-elliptic Weierstrass functions (291. Without the constant in v, these are the 
usual (non-normalized) Bloch frmctions; nevertheless, we choose here a different normalization for both 
of them, i.e., u without the usual constants and v with a prefactor that depends on <5. This is so for having 
the minimum necessary factors in order to fulfill the orthogonality condition W(m, v) = 1, which can be 
shown using Eqs. (SOI and dSTIl from the Appendix. 

Note thatis defined by the relation m(x + T) = ^u[x) and then= exp[2^^((u) —2a;^(^)]. Besides, 
by expressing it as /3 = e”^, then k = 2i[co[,{S) — <5f(a;)] (up to an additive multiple of 2ni) which is 
known as the quasi-momentum 13^ . From the analytical form of this quantity, we can obtain the allowed 
and forbidden bands by calculating if k takes real or complex values, respectively. The displacement S 
and the factorization energy e are related by Il27l : 

e + - p{S). (72) 

It is worth pointing out that we are using one Bloch state to perform the SUSY transformation, even 
when these states carmot be used for non-confluent SUSY transformations (a linear combination of the 
two Bloch solutions must be used in that case Il3l f35l ). Nevertheless, one of the advantages of the conflu¬ 
ent algorithm is that it can accept a larger set of eigenstates to perform the transformation, for example a 


(70) 

(71) 
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Figure 3: Some SUSY partner potentials and their added states obtained from the single-gap Lame potential 
using the 3-confluent SUSY transformations. The plots correspond to the original potential Vg (*) with m = 0.5 
(blue) and two modified ones: V 3 (x, 0.5,—0.2,0.1,0.01,0.01) (orange) and V 3 (x, 0.5,—0.5,0.1,10,0.01) (purple) 
with the format 'V^{x, m, e, Cj, Dj, D 2 ). As m = 0.5 and e = {—0.2, —0.5}, the transformation energy is in the 
lower forbidden band. (Color online). 


state that diverges in one side of the domain. These functions cannot be used in the non-confluent cases 
and this is exactly the situation for the Lame potential. 

We are going to calculate next its parametric derivative with respect to e, for which we will employ 
some identities for the elliptic functions shown in the Appendix (see Eq. (l80t ). Then, using the chain 
rule and Eq. (l72l we get: 

deU = (deS) dsU = -{dsp)~^dsU. (73) 

An explicit calculation produces: 

dgu = [^{x -|- (5 -|- co’) — ^{3 co') xp(^S)]u, (74) 

With these equations we are now able to calculate the Wronskian using symbolic computation soft¬ 
ware, either with Eqs. ll29ll or Eqs. ll30ll . In this case, Eq. (l30t is easier as it is written in terms of simpler 
Wronskians. The result is too long to be written here explicitly; however, we plot several examples of 
the potentials and their added states in Fig. |3l We obtain a 4-parametric family of potential with this al¬ 
gorithm, i.e., e, Cl, Di, D 2 . As in the previous case, only a subset of these solutions will be non-singular, 
unfortunately in this case it is not easy to determine a priori all the constraints for the four parameters, 
although we know that e i Sp(H). 

The calculations for the 4-confluent case are much more complicated, but thanks to the fact that this 
method depends only on calculations of derivatives with respect to x and e, symbolic software is still 
able to calculate analytically the SUSY partner potentials and the added states for the Lame potentials. 
Using Eqs. (l37l l, l l38ll and (|39] | we calculated the new potentials and with Eq. (l40l their physical added 
states. The results are again too long to be written here but we plot some results in Figs.|4]and0 Given 
that we are using the confluent algorithm, we can use transformation energy in any of the forbidden 
bands, including the band gap, which can only be achieved by traditional SUSY transformation by 
a limiting procedure and leads to the same family of potentials as the confluent case (see Fig. 4 in 
Ref. 13 and Fig. 8 in Ref. l35ll ). We present a a transformed SUSY potential and their added states 
using transformation functions from the band gap and from the lower forbidden band in Figs. |4] and |5l 
respectively. 

We must remark that even though we are only plotting V 3 , ip^, V 4 , ip^ for the Lame potential, we also 
have exact analytical expressions for all of them, they are just too long to be written here. 
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Figure 4: Some SUSY partner potentials and their added states using 4-confluent SUSY transformations de¬ 
parting from the single-gap Lame potential. The plots correspond to the original potential Vg (^) with m = 0.5 
(blue) and two modified ones; V 4 (x, 1.25,0,0,1,0,0,0) (orange) and V 4 (x, 1.45,0,0,0,0,1,0) (purple) with the 
format Va (x, m, e, Ci, Ct, Di, Do, Dg). In this case the transformation energy belongs to the band gap as m = 0.5 
and e = {1.25,1.45}. (Color online). 



Figure 5: Some SUSY partner potentials and their added states using 4-confluent SUSY transformations de¬ 
parting from the single-gap Lame potential. The plots correspond to the original potential Vg (^) with m = 0.5 
(blue) and two modified ones: V 4 (x, 0.5,—1,1,—2,—1,2,3) (orange) and 14 (x, 0.5,—1,0,0,0,0,—1) (purple) 
with the format V 4 (x, m, e, Cg, C 2 , Di, D 2 , D 3 ). The transformation energy is in the lower forbidden band as 
m = 0.5 and e = — 1. (Color online). 
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7 Conclusions 

In this work we have developed the general framework for the Ic-confluent SUSY transformation using 
a Wronskian formula which involves only derivatives with respect to x and to the factorization energy 
e instead of integrals, as in the usual fc-confluent case. In order to do this, we have explicitly worked out 
the second-, third- and fourth-order cases. Then we have proven the general formulas for the fcth-order 
case by mathematical induction in Sec. |5] 

For the 3- and 4-confluent cases we have also developed some simplifications of the respective Wron¬ 
skian formulas. We believe that these formulas (see Eqs. ll29ll - (l30ll and (l36l - (l38t l can be very useful 
because they reduce the number of calculations necessary to obtain the transformed potential. 

Then, we have applied this formalism to two potentials. The first one is the free particle potential, 
for which the 3-confluent transformation has been calculated before ESI; our calculations here work as 
a check for our algorithm. We foimd the expected results already found in the literature. Our formula 
gives us a 4-parametric family of potentials. Furthermore, we calculated the 4-confluent SUSY partner, 
which gives us more freedom to deform this potential initially through six parameters, although after we 
find the conditions to restrict to the non-singular subset, we arrive again to the Poschl-Teller potential, 
exactly the same results as in the 3-confluent case. The reason is that any non-singular potential obtained 
from the free particle with only one added level will always have this form ll4ll5l l23ll2H . 

The second one is the single-gap Lame potential, for which the calculations become much more 
difficult due to the natural complex nature of the elliptic functions such that most of them are complex 
and only for special values of their parameters they become real (up to numerical errors in the imaginary 
part of the order of 10”^®). Nevertheless, we were able to show the possibilities of this method, as for 
this potential we were able to calculate its 3- and 4-confluent SUSY partners potentials analytically. 

We believe this method can be applied to more potentials, as several quantum systems can be solved 
in terms of the confluent hypergeometric function i Fj, and there is a work where their parametric deriva¬ 
tives have been worked out l39l . 

We must remark the strengths of this method. First, it is a confluent SUSY transformation, which 
means that it can be applied to a more general class of transformation function than the usual SUSY or 
Darboux transformation, i.e., using some non-physical states (for the Lame equation we were able to use 
Bloch states). Second, it is the only way to add only one level above the ground state energy. This case 
can be obtained as a confluent limit from the usual second-order SUSY transformation. Third, because 
we have developed the higher-order algorithm and therefore we have now further tools for deforming 
the potential, because the higher the fc-confluent transformation, the higher the number of parameters 
that the new potentials will depend on. Also, we perform this higher-order transformation in a single 
step. Fourth, the calculation is done through parametric derivatives and Wronskians, unlike the usual 
confluent case that uses integrals, that most of the times are not easy to solve (they are even unsolvable). 

Furthermore, in the case of the Lame potential the possibility of adding a new level in the band gap 
is significant, because this is the simplest way to model impurities in solid-state physics. The transitions 
between the valence and the conduction bands are easier when a semi-conductor has impurities, this 
can be modeled by adding an intermediate physical level in the band gap Il38l . An example of this new 
bound state is in the plots of Fig.|4]and the spectrum will look like the one in Fig. [6] 

Finally, there is still perspective work related to the non-singular subfamilies of potentials. Also, it 
seems likely that there is an the underlying structure for the Wronskian formula for Uconfluent SUSY 
which is worth of further study. 
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Figure 6: Spectrum of the deformed single-gap Lame potential. The usual structure of the single-gap Lame 
potential has one more bound state (the white line in the middle of the band gap), in solid-state physics it is 
called localized impurity state. 


Appendix 

By using the following Wronskian identity 

Wif,hf) = {d,h)f, (75) 

which is valid for two differentiable but otherwise arbitrary functions / and h, it is straightforward to 
show that the Wronskian formula for the SUSY transformation given by Eq. iQI is preserved for the 
confluent case. 

Using that u and v fulfill the Schrodinger equations (IT2ll and dT^ and from the Jordan block of length 
3 in Eqs. (l23ll after some algebra we get 

W(mi,M2 ,M3) =UiW{ui,U3)-U2W{ui,U2). (76) 

Similarly, for the Jordan block of length 4 in Eqs. ll33ll we get 

W(mi,M2, 1/3,14) = W(mi,M2)[W(mi,M4) -L W(m 2,M3)] “ W^Ju^Ms). ( 77 ) 

Deriving W{u,v) = 1 with respect to e we obtain 

W{u,deV) = W{v,deU), (78) 

and deriving again 

W(deU,deV) = i[W(M,deU) - W(v,d^u)j. (79) 

Also, there are some relationships between <7(x), ^(x), and p(x) and its derivatives that can be found 
in Ref. 1291 and are shown next: 

dxcr(x) = (7(x)((x), (80a) 

d:c^(x) = -p(x), (80b) 

(80c) 

Erom Ref. l28ll we also have 

+ (81a) 

<7(zi -h Z2)<7 (Zi - Z 2 ) = -cr(zi)£7-(z2)[p(zi) -e(Z 2 )]. (81b) 

and the fact that p is an even function and u and ( are odd. 
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